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We develop a quantitative description of giant asymmetry in reflectance, recently observed in 
semicontinuous metal films. The developed scaling-theory based technique reproduces the spectral 
properties of semicontinuous composites, as well as provides insight into the origin of experimentally 
observed loss, reflectance, and transmittance anomalies in the vicinity of the percolation threshold. 
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I. INTRODUCTION 

Research into the optics of semicontinuous metal- 
dielectric films has been enjoying sustained interest 
due to a unique combination of novel physics and the 
practical applications offered by such composites. It 
has been demonstrated, both theoretically and exper- 
imentally, that the electromagnetic (EM) response of 
these structures is dominated by a non-trivial inter- 
play between Anderson-localized and delocalized sur- 
face plasmon polaritonsji"— This results in unusual op- 
tical properties that include: greatly enhanced absorp- 
tion, giant intensity fluctuations of local EM fields, gi- 
ant local chiral response, and strongly enhanced optical 
nonlinearities.— In a related context, it has been re- 
cently shown that the reflectance of semicontinuous sil- 
ver nanocomposites, chemically deposited on glass sub- 
strates, strongly depends on the direction of incident 
lightjii In particular, the reflectance of such a system 
irradiated from the substrate/film interface side can dif- 
fer by as much as 15% from its reflectance given film/air 
side incidence (see FiglT|). Moreover, this large asymme- 
try in reflectance has been found to be extremely broad- 
band, spanning most of the visible frequency spectrum. 
For comparison, the reflectance asymmetry of thin, con- 
tinuous silver films does not exceed 3% when measured 
over the same range of optical frequencies, and does not 
exhibit any broadband characteristics. It has been sug- 
gested that the origin of this large broadband asymmetry 
is in the enhanced optical absorbance which is often seen 
in percolation-type systems. Here we develop a quanti- 
tative description of the observed phenomenon. 

The geometry of the system described in this 
manuscript is shown in Figure [TJ We approximate the 
silver percolation film as a uniform material with thick- 
ness d [in our calculations d = 50nm (see section III)]. 
The microstructure of the film is characterized by the 
surface metal filling fraction p ranging from p = for 
a bare glass substrate to p = 1 for a substrate which is 
fully covered with metal. At the critical value p — Pc, 
known as the percolation threshold, the do conductivity 
response of the entire random metal-dielectric compos- 
ite undergoes an insulator-conductor phase transition. 
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FIG. 1: (Color online) General layered structure composed 
of a silver percolation film clad by air to the left and glass 
to the right. Incident light may come from either the air or 
substrate side as shown. Both air and glass regions are taken 
to be semi-infinite. 



The unique optical properties of our films are manifest 
particularly in the vicinity of the percolation threshold, 
and we therefore adopt the conventional description of 
the response as function of the parameter p — pc for the 
purpose of both modeling and data analysis. All theoret- 
ical and experimental results in this work use pc = 0.6, 
which correlates well with two dimensional site percola- 
tion on a square lattice {pc — 0.593) 

As mentioned above, the reflectance Ri of the com- 
posite film, measured using light impinging from the 
air/metal interface, strongly differs from i?2 - the re- 
flectance measured when light is incident from the sub- 
strate/film side. Since the transmittance of our system, 
as the transmittance of any non-chiral homogeneous film 
is symmetric {i.e. Ti = T'2)^'^'^*, the asymmetry in re- 
flectance Ai? = i?i — i?2 directly reflects the asymmetry 
in losses As we show below, in contrast to vacuum- 
deposited percolation films, (i) Ai? as well as the com- 
puted combined losses exhibit a local minimum at p ~ pc, 
(ii) Ai? exhibits broadband response in the vicinity of 
p — Pc — ±0.05 , (iii) the reflectance exhibits a local max- 
imum in the vicinity of p ~ Pc, and (iv) the transmittance 
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exhibits a local minimum near p c:i pc- 



II. PERCOLATION FILM SYNTHESIS AND 
CHARACTERIZATION 

Semi-continuous silver films with controllable filling 
fractions were deposited on microscope slides using a 
modified ToUen's reaction as described previouslyj ^^i^^ 
The amount of silver deposited on the substrates was 
controlled by monitoring deposition times, with reactions 
ranging between 1-6 h. Ensuing deposition, coated sub- 
strates were rinsed with ultrapure water before being 
dried with filtered air and stored under nitrogen until 
tested. Figure [5] shows a scanning electron micrograph 
of a typical film with filling fraction p ~ 0.52. These 
chemically deposited films appear as highly disordered 
polycrystalline aggregates, with large grain size distribu- 
tions. In addition, we note the non-uniform coating of 
the substrates by the metal, resulting in highly discon- 
tinuous morphologies. 




FIG. 2: (Color online) Scanning electron micrograph of a 
chemically deposited silver film with metal filling fraction 
p ~ 0.52. The scale bar is 500nm. 

Optical refiectance and transmittance spectra were col- 
lected using a spectroscopic optical microscopy setupi^i 
A tungsten-halogen white-light source was used to illumi- 
nate the samples through an inverted microscope whose 
output was imaged on the entrance slit of a _F = 320mni 
spectrometer with a resolution of 0.5 nm. A lOx (0.25 
N.A.) objective was used to collect the normally incident 
light for spectroscopic imaging onto a liquid nitrogen- 
cooled detector. A high-rcflcctance mirror (Newport 
Broadband SuperMirror, R > 99.9%) was used to nor- 
malize all signals. In addition, data from ~ lOOO/im 
across the films were averaged to obtain each final trace, 
in order to eliminate spurious effects from local inhomo- 
geneities in the rough films. The spectral response of the 



film shown in Figl2]is depicted in Fig|3l 



III. REFLECTION, TRANSMISSION, AND 
ABSORPTION OF RANDOM PERCOLATION 
COMPOSITES 

Many metal-dielectric composite systems are de- 
scribed by conventional effective medium techniques 
(EMTs)i2r— by representing the composite as an effec- 
tive homogeneous layer which successfully models the 
system's average optical properties. However, it is known 
that the optical properties of these films close to the 
percolation threshold cannot be adequately described by 
EMTs,2ii2^ The reason for the consistent failure of EMTs 
in this case is two-fold. First, although the dimensions 
of the components in percolation films are much smaller 
than the free-space wavelength, the optical properties of 
the composites are dominated by the dynamics of res- 
onant clusters that can be comparable in size to the 
wavelength. Second, as result of a dc metal-dielectric 
phase transition, the effective parameters of the percola- 
tion films in vicinity of Pc become scale-dependent and 
therefore cannot be described by quasi-static effective 
medium models. Although some percolation films have 
been successfully described in terms of generalized Ohm's 
law (G0L)2^"— , straightforward extensions of GOL for- 
malism to our system are not consistent with our exper- 
imental observations. 
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FIG. 3: (Color online) Measured reflectance (red diamonds), 
transmittance (blue boxes), and absorbance (green circles) 
as function of incident wavelength for measured metal filling 
fraction p ~ 0.52. Solid lines represent the results of scaling 
theory calculations. 

The only technique that explicitly accounts for the dc 
conductivity phase transition of percolation systems is 
known as scaling theory i '^^^^^^^^ — In this technique, the 
conductivity of the film is assumed to be explicitly de- 
pendent on the size of the cluster L over which it is mea- 
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sured. Specifically, the average conductivity of a conduc- five cluster of size L is given by 
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and that of a dielectric (insulating) cluster is 
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The expressions above explicitly assume that the con- 
ductivity of the conductive component of the film is given 
by the Drude model, 
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where Odc is the dc conductivity, t is the electron relax- 
ation time, and w is the angular frequency of the incident 
light. The ac response of a dielectric film component is 
equivalent to that of a capacitor, 
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where Cq is the average capacitance between neigh- 
boring metal clusters. The parameters (Tdc,T, Cq 
and Ci . . . C4 coefficients are uniquely determined by 
the composition and micro-geometry of the percolation 
film. The critical exponents for 2D percolating films are 
^ — s — l.'i and V — 4/3 ,21,22 ,29 p <^ percola- 
tion films are governed by dielectric conductivity, which 
is dominated by the capacitance coefficients C3 and C4. 
For p :s> Pc metallic conductivity (governed by Ci, C2) 
dominates the optical properties of the system. 

Despite the scale-dependence on the microscopic and 
mesoscopic levels, the percolation film appears homoge- 
neous when measured over a significantly large area. The 
transition from the scale-dependent to the homogeneous 
dc response occurs at the scale known as the correlation 
length, ^, that characterizes the typical cluster size. In 
the vicinity of the percolation threshold, the correlation 
length diverges as 



P-Pc 
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(5) 



The constant ^0 represents the smallest metal cluster size, 
which occurs at p — ^ 0. 

At finite frequencies, the oscillatory motion of elec- 
trons within conducting clusters leads to the length scale 
correction of the homogeneous response of the system. 
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where 9 = 0.79, Bq = 4.0, and the free space wavelength 
is given by Aq.^^'^^ 

The ac conductivity of the percolation films, calculated 
using the expressions above can be directly related to an 
effective film index, which along with the film thickness 
d can be used to determine the macroscopic optical prop- 
erties of the film, including R and T. In our calculations, 
we use the technique introduced in Rcf.[21]. In this ap- 
proach, the optical properties of the film are calculated 
as a weighted average of conductive (dielectric) film con- 
tributions, where the average conductivities are given by 
Eqs.([T]) and ^ respectively to yield 
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[fR^Az^rn) + (1 - f)R,,a{z<Jd)]P{z}dz (8) 



where the parameter 
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is the metal occupation probability, 
metal concentrations, p < p, 
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pation probability / — >■ 0. When p > pc 

the occupation probability / — >■ 1. For intermediate sur- 
face metal concentrations centered at p = Pc with full- 

width Ap — 2pc ij^j , the occupation probability 

varies linearly as a function of p from unoccupied (/ = 0) 
to occupied (/ — 1). As shown by the above inequali- 
ties, the range of metal surface coverage values for which 
scaled metal and dielectric optical properties are aver- 
aged depends non-trivially on the correlation length, ap- 
plied frequency, and film geometry. The function P{z) 
gives the distribution of the conductivities of conductive 
[dielectric] clusters around their mean values given by 
Eq.© [Eq.©]. Following Ref.[2l| and we assume 
that P{z) is adequately described by a log-normal dis- 
tribution function with standard deviation of asd = 0.3. 
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Integrating over all scaled conductivities averages out the 
length dependent optical conductivity and allows for per- 
colation films to be modeled by the contributions from 
planar homogeneous constituent layers. 

The homogeneous-layer optical properties are given 
by^ 
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where the glass substrate index is ng = 1.5166, the ef- 
fective film index is Uf = ^1 + Airier /uj, and the phase 
parameter is 4> = exp(i^nf(i). 

As noted, all critical exponents in the expressions 
above are universal for all 2D percolating networks, while 
the parameters Cq . . . C4, Udc, t, and ^0 are unique for a 
given percolation filmi^ii^ In our calculations, we use 
tJdc = 2.574 X 10^''sec~^, frequency dependent relaxation 
time^ l/r = I/tq + Poj^, where tq = 3.0fs, /3 = 0.2fs, 
Co = 0.5, Ci = C2 = 0.046, C3 = 0.028, C4 = 0.055, and 
^0 = 2nm. 



IV. COMPARISON WITH EXPERIMENTAL 
RESULTS AND DISCUSSION 

A comparison of the experimentally obtained spectral 
response of the silver films with the predictions of scaling 
theory is shown in FiglH It is seen that both the broad- 
band nature of the reflectance asymmetry and its non- 
monotonic behavior near the percolation threshold are 
are well reproduced by the theoretical model, as demon- 
strated in FigUl We note however, that the model fails 
for the large metal concentrations p — > 1— li^ where the 
structure of the composite becomes substantially three- 
dimensional and cannot be treated as a thin homoge- 
neous film. To further illustrate the robustness of the 
presented technique we show in Fig|S] a comparison of 
experimentally measured values of Ri and Ti, as well as 
the losses (computed as ^1 = 1 — Ri — Ti,) with our 
theoretical model. As mentioned before, both theoreti- 
cal and experimental results clearly show that despite a 
strong reflectance asymmetry, the transmittance of the 
films remains symmetric. Therefore the asymmetry in 
reflectance is directly related to the asymmetry in losses 
{AR = i?i - i?2 = A2 - Ai)M 

We now examine the loss (or, alternately the re- 
flectance) in more detail. We note that several previous 
experiments as well as theoretical models have observed 
absorption maxima in the vicinity oi p = Pc. In contrast, 
our experimental data clearly demonstrate a local mini- 
mum in losses near the percolation threshold, as seen in 
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FIG. 4: (Color online) Points represent the measured change 
in reflectance (Ai? — Ri ~ R2) for various incident wave- 
lengths. Black circles 500nm, green triangles 600nm, red 
boxes 700nm. Corresponding colored solid lines (black solid 
500nm, green short-dashed 600nm, red long-dashed 700nm) 
represent the results of scaling theory reflectance calculations. 
The inset shows the change in reflectance over the entire sur- 
face coverage range. Note that the 2D scaling model fails 
for large metal concentrations, where the three-dimensional 
structure of the composite dominates the optical response. 



FiglSJe). Additionally we observe a local maximum in 
reflectance and a local minimum in transmittance near 
p = Pc, as seen in Figl5l[a,c). 

We suggest that this anomalous behavior stems from 
the dramatically reduced correlation length in our 
solution-derived percolation films. The smallest metallic 
particle size produced in our experiments is on the order 
of 2 nm, in contrast to 10 nm reported in Ref . [22j . In 
addition, the optical response of solution-derived metals 
is typically aff'ected by the reduced electron mean free 
path'^^, which may further reduce the effective particle 
size, described by the parameter ^o- 

Figure [5] demonstrates the evolution of optical proper- 
ties of percolation systems when the correlation length is 
reduced, corresponding to a change in ^0 from 10 nm to 
2 nm. It is clearly seen that at ^0 — 2nm the absorp- 
tion reaches a local minimum in the vicinity of p = Pc, 
while at ^0 = lOnm, the system recovers the absorption 
maximum at p ^ Pc, as observed in previous studies. 

The dramatic difference between correlation lengths 
in our solution-derived filmsii^i^ and other vacuum- 
deposited counterparts2ii2^i24 is consistent with the dif- 
ference in fabrication technique: while thermal deposi- 
tion under vacuum typically yields uniform metal films 
with almost perfect Ohmic contacts between adjacent 
grains, solution-based deposition is routinely associated 
with quantitatively weaker contacts between conducting 
grains. The latter results in films with reduced electron 
mean-free-paths and lower correlation lengths. 
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FIG. 5: (Color online) Points represent the measured (a) re- 
flectance, (c) transmittance, and (e) loss from the air side as 
a function surface coverage fraction. Connecting lines are a 
guide for the eye. Calculated (b) reflectance, (d) transmit- 
tance, and (f) loss when the correlation length parameter is 
^0 = 2nm (solid line), ~ 5nm (dashed line), and = lOnm 
(dotted line). For all graphs the incident wavelength is 700nm. 



V. CONCLUSION 



We have developed an analytical description of the 
phenomenon of broadband asymmetric reflection in per- 
colation composites. The developed technique, based on 
scaling theory, is not only capable of describing the spec- 
tral response of our films, but also explains that the re- 
duced correlation length in our solution-derived compos- 
ites is the primary cause of the experimentally observed 
anomalous optical properties near the percolation thresh- 
old. Our work demonstrates that the correlation length 
is an important factor that fundamentally affects the op- 
tical properties of percolation composites in the vicinity 
of the percolation threshold. 
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